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ABSTRACT 

We present a two-dimensional mapping of the gas flux distributions, as well as of the 
gas and stellar kinematics in the inner 220 pc of the Seyfert galaxy NGC 2110, using K-band 
integral field spectroscopy obtained with the Gemini NIFS at a spatial resolution of ~ 24 pc 
and spectral resolution of « 40 km s -1 . The H2 A 2.1218/mi emission extends over the whole 
field-of-view and is attributed to heating by X-rays from the AGN and/or by shocks, while 
the Bry emission is restricted to a bi-polar region extending along the South-East-North- 
West direction. The masses of the warm molecular gas and of the ionized gas are Mr 2 ~ 
1.4 x 10 3 M© and Mr ii « 1.8 x 10 6 M©, respectively. The stellar kinematics present velocity 
dispersions reaching 250kms _1 , and a rotation pattern reaching an amplitude of 200kms _1 . 
The gas velocity fields present a similar rotation pattern but also additional components that 
we attribute to inflows and outflows most clearly observed in the molecular gas emission. 
The inflows are observed beyond the inner 70 pc and are associated to a spiral arm seen in 
blueshift to the North-East and another in redshift to the South-West. We have estimated a 
mass inflow rate in warm molecular gas of « 4.6 x 10“ 4 M©yr -1 . Within the inner 70 pc, 
another kinematic component is observed in the H 2 emission that can be interpreted as due to 
a bipolar nuclear outflow oriented along the East-West direction, with a mass-outflow rate of 
« 4.3 x 10 -4 M© yr -1 in warm H 2 . 

Keywords: galaxies: individual (NGC 2110) - galaxies: Seyfert - galaxies: ISM - infrared: 
galaxies - galaxies: kinematics and dynamics. 

clear activity in galaxies suggests that the gas might be stream¬ 
ing along these structures toward the central region to feed 
the SMBH inside ([Ferrarese & Fordll2006l ; ISomerville et al]l2008l : 
Kor mendv & H0II2OI 3h . The dust mass at kp c scales in AGNs is 
estimated to be in the range 10 5 - 10 7 M o dSimdes et Idl 120071: 
Mart ini et alJl2013h and the associated large amounts of molecu- 
lar gas (10 7 to 10 9 M 0 ) reinforces the importance of searching for 
signatures of molecular-gas inflows within the inner kiloparsec. 

On the other hand, r ecent studies have also found massive out¬ 
flows of molecular gas dSakamoto et al.1 l20ld : lAalto et ahl 120121 : 
IVeilleux et alJEoIA in particular, in Luminous Infrared Galaxies 
(LIRGS) and Ultra Luminous Infrared Galaxies (ULIRGS). 

Observations of nearby active galaxies using spatially re¬ 
solved IFS provide constraints for the feeding and feedack pro¬ 
cesses. The study of these processes in nearby AGNs have been 
carried out using IFS at near-IR bands, which contain emission 
lines from molecular (tracer the AGN feeding) and ionized (tracer 
of the AGN feedback) gas. Most of the observations are performed 
with adaptive optics using 8-10 m telescopes at spatial resolutions 
ranging from a few parsecs to tens of parsecs. The main results 
of these studies are: (i) In general, the molecular and ionized gas 
present distinct flux distributions and kinematics. The H2 emission 


1 INTRODUCTION 

A necessary condition for the onset of nuclear activity in galax¬ 
ies is the presence of gas to feed the nuclear Supermassive Black 
Hole (hereafter SMBH) iHeckman & Besdl2014h . Evidence for the 
presence of large quantities of gas in the inner few hundred par¬ 
secs of nearby active gal axies has been indeed found in a number 
of imaging studies (e.g. iMalkan et al.1 1998). These studies reveal 
the presence of dusty structures, in the jfor m of compact nuclear 
disks, filaments and spirals. In ISimoes et al.1 i2007lk it was found 
that, in the particular case of an early-type sample of active and 
non-active galaxies, all the active galaxies had dusty structures in 
the inner few hundred parsecs, while only 25% of the non-actives 
had such structures, implying that the dust was a tracer of the gas 
feeding the active nucleus. The excess of dust in the active galaxies 
relative to the non-active ones was later confirmed by a study us¬ 
ing Spitzer observations to calculate the dust mass of these regions 
jMartini et al.112013l) . As the inferred gas mass around the active 
galactic nuclei (AGN) are much larger than necessary to feed the 
SMBH inside, these structures are reservoirs of material where new 
star formation may be triggered. 

This connection between the presence of dusty spirals and nu- 



























2 Diniz, Riffel, Storchi-Bergmann & Winge 


is dominated by rotational motion and shows lower velocity dis¬ 
persion values than the ionized gas. From the gas excitation and 
kinematics, we have found that the emitting H 2 gas (hot molecu¬ 
lar gas , as it originates in gas at « 2000 K. IStorchi-Bergmann et al.l 
(12009b ) is located in the plane of the galaxy, presenting a ve¬ 
locity field dominated by rotation. In some cases, in the inner 

« _ 30 pc of the galaxy, compact molecular disks dHicks et al] 

201 3l: Mazzalav et all 1201 4 IRiffel & Storchi-Bergmannl 1201 lal : 


Scho nell et alJ 12013h a nd gas inflows along nuclear spirals (e.g. 
Riffel et al.ll2008l . 120131) are found. The inflow rates in cold molec¬ 


ular gas range from a few tenths to a few solar masses per year 
(Riff el et al. 2008L 20091: Davies et al. 20091: ISanchez et all [2009: 


IRiffel & Storchi-Bergmannl201 lal :iRiffel et al.1201317 . (ii) Although 

the ionized gas kinematics frequently presents a rotating disk com¬ 
ponent similar to that seen in H 2 for most galaxies, it shows a 
large contribution from an outflowing component for some objects 
(Riffe l et akll2006[ 20091 120131 1201 4bl: Riffel & Storchi-Bergman 


201 lalld: LStorchi-Bergmann et al.l l2010l: ISchonell et al 1201 


i 


Barbosa et al. 


201 4) . These outflows are most easily seen in the 


[Fe 11] 1.257 /d m and 1.644yum emission lines, present in the near-IR 
spectra. In IBarbosa et al.l d2014h it is shown that these lines seem 
to be better tracers of AGN outflows than the previously used [O 
hi] A 5007 A emission line. This is due to the fact that the [Fe 11] 
emission is more extended than that of the [Om], reaching out 
to the partially ionized region (beyond the fully ionized region). 
Mass-outflow rates in ionized gas ha ve been estimated to be inthe 
range from 0.05 to AM^yr -1 ( Riffel & Storchi-Bergmannl 2011c ; 
Storc hi-Bergmann et al.l 120101 : IRiffel et al. 20131 : IBarbosa et al] 


20141) . 

In this work we present adaptive optics assisted K-band IFS 
of the inner 3" x 3" of the galaxy NGC2110 obtained with Gem¬ 
ini NIFS. The results shown here are part of a larger project of the 
AGNIFS (AGN Integral Field Spectroscopy) group to map inflows 
and outflows in a complete sample of AGN, selected by their X- 
ray emission. NGC2110 is an early type SO galaxy which hosts 
a Seyfert 2 nucleus, located at a distance of 30.2 Mpc, assuming 
a redshift z = 0.007789 (NEE0), where 1" corresponds to 146 pc 
at the galaxy. Due to its proximity, it has been already observed 
in several wavebands, from radio to X-rays. NGC2110 is one of 
the brightest nearby Seyfert galaxies in hard X-ray s (2- 10keV) 
and is classified as a narrow-line X-ray galaxy jBradt et al.lll978h . 
Its luminosity is comparable to that observed for Seyfert 1 nuclei 
and soft X-ray emission was observed extending up to 4" to the 
north of the nucleus (We aver et al.ll 1995k Dust lanes and asymme¬ 
tries a re seen up to 4" west of the nuc leus in optical jMalkan et al.l 
1 1998b and near-IR dOuillen et all 1999b bro ad band images obtained 
with the Hubble Space Telescope (HST). IStorchi-Bergmann et akl 
( 1999 ) obtained long-slit near-IR spectroscopy of the inner ~ 10"of 
NGC2110 along the major and minor axes of NGC2110 finding 
strong emission of both [Fen] and H 2 , as well as signatures of dis¬ 
tinct kinematics for these two species. Excitation by X-rays was 
favored for H 2 , while shocks were concluded to be also important 
for excitation of the [Fe n]. 

Observations in radio reveal extended emission in a S-shaped 
structure, with an extent of « 4" along the N orth-South direction 
JlJlvestad & Wilsonlll983l : iNagar et al. Il999b. Circum nuclear ex¬ 
tended gas emission was studied bv Evans et al.ld2006b using Chan¬ 
dra, HST and VLA observat ions. Long-slit spectroscopic observa¬ 
tions of the inner 10" jGonzalez Delgado et al.ll2002 ; Ferruit et al.l 
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12004 have shown that the gas velocity field is asymmetric relative 
to the nucleus at distances larger than 1", while within the inner 
1" a blueshift excess is observed, suggesting the presence of a nu¬ 
clear outflow. [Rosario etal . ( 2010 ) have shown that the outflow is 
oriented at a position angle (PA) offset by « 40° from the PA of the 
radio jet and suggested that the gas is ionized by the central AGN 
and not by shocks. 


The gas kinemati cs was first investi gate by_ Wilson & Baldwinl 
dl985ab and I Wilson. Baldwin & Ulvestadl (1 1985b . who reported a 
rotation pattern in the ionized gas kinematics, similar to those ob¬ 
served for other spiral galaxies, suggesting that the gas is in ro¬ 
tation in the galaxy disk. They found that the rotation center was 
located 1.7" south of the photometric nucleus and observed asym¬ 
metric profiles in the [0111] t o the southeas t. In the long -slit near-IR 
spectroscopic observations of IStorchi-Bergmann et al.l dl999b . they 
also found a displacement of the kinematic nucleus, but smaller 
than the optical, thus attribuited at least partially to obscuration, 
which is smaller in the near-IR than in the optical. Recent op¬ 
tic al IFS with the Gem i ni Mu lti-Object Spectrographs (GMOS) 
by ISchnorr-Muller et all d2014b revealed the presence of four dis¬ 
tinct kinematic components in the emitting gas: a cold gas disk, 
a warm gas disk, a compact nuclear outflow and a high latitude 
cloud identified as “the northern cloud”. They also observed excess 
blueshifts and redshifts in the cold gas disk kinematics relative to 
rotational motion that were attributed to inflows and outflows.They 
have shown that the apparent asymmetry of the gas velocity field is 
due to the superposition of the contribution of the different compo¬ 
nents; when this is taken into account, there is no more asymmetry 
in the velocity fields. 


Studie s of the nuclear stellar population by 
Gonzalez Delgado e t al ] teomli using optical continuum im¬ 
ages did not show signatures of a young stel lar population (age < 
1 Gyr). On the other hand, iDurre & Mould I d2014h used near-IR 
(Z, J and H bands) IFS of the inner 07 8 obtained with the OSIRIS 
instrument at the Keck Telescope at an angular resolution of 70 
mas, and reported the detection of four star forming clusters at 
distances smaller than 073 from the galaxy nucleus. They also 
found that the ionized gas (in particular [Fen]) surrounding 
the star clusters is being excited by strong outflows associated 
with recent star formation and derived star formation rates of 
0.1 2-0,26 M 0 yr -1 , The new fin dings by IDurre & Mould I d2014b 
and ISchnorr-Muller et all d2014b . relative to previous studies can 
be attributed to the better spatial coverage of IFS to study the 
central region of galaxies when compared to long-slit studies. 
iDurre & Mould 1 ( 2014 ) also calculated the SMBH mass using 
the He I velocity profile considering a simple Keplerian rotation 
in a region with radius 56 pc and found a mass for the SMBH of 
« 4 x 10 8 M o . Using the M-cr relationship they find a value of 
~ 3 x 10 8 M o while with the iGraham et al . ( 1201 lb updated rela¬ 
tionship for elliptical galaxies they derive a value of * 5 x 10 7 M o ; 
a lower value for the SMBH mass was also obtained using the 
[Fen]velocity field, probably due to the presence of outflows from 
the clusters. 


This paper is organized as follows. In Sec. [2] we describe the 
observations and data reduction procedures. The results are pre¬ 
sented in Sec.[3]and discussed in Sec. [4] We present our conclusions 
in Sec. [5] 
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2 OBSERVATIONS AND DATA REDUCTION 

The near infrared observations of NGC 2110 were obtained with 
the integral field spectrograph NIFS (iMcGregor et al.1 l2003h on 
the Gemini North telescope in two distinct nights, under the pro¬ 
gramme GN-2010B-Q-25. In the first night, five exposures of 600 
seconds were obtained for the galaxy, 3 observations with the same 
exposure time for the sky and 3 exposures of 18 seconds each for 
the telluric star. In the second night, only one exposure of 600 sec¬ 
onds for the galaxy and for the sky and 4 observations of 15 seconds 
for the telluric standard star were obtained. 

The data reduction was performed using specific tasks con¬ 
tained in the nifs package which is part of gemini iraf package. 
The reduction process followed standard procedures such as trim¬ 
ming of images, flat-fielding, sky subtraction, wavelength and spa¬ 
tial distortion calibrations. The telluric absorptions were removed 
using the spectrum of the telluric standard star and the flux cali¬ 
bration was performed by interpolation a black body function to 
the spectrum of the telluric standard star. We combined the individ¬ 
ual datacubes into a single datacube using the gemcombine task of 
the gemini iraf package and used the sigclip algorithm to eliminate 
the remaining cosmic rays and bad pixels, using the peak of the 
continuum as reference for astrometry for the distinct cubes. The 
final data cube covers the inner 3"x3"(~440x440pc 2 ) region and 
contains ~ 3900 spectra, each spectrum corresponding to an angu¬ 
lar coverage of 0705 x0705, which translates into ~7x7pc 2 at the 
galaxy. 

The spectral resolution of the data is about 3.4A at 2.3 pm, 
measured from the full width at half maximum (FWHM) of the arc 
lamp lines and corresponds to a velocity resolution of 44kms _1 . 
The angular resolution is 0715, corresponding to « 22 pc at the 
galaxy, and was obtained from the FWHM of a Gaussian curve fit¬ 
ted to the flux distribution of the standard star. A similar FWHM 
is found by fitting the unresolved nuclear K-band continuum emis¬ 
sion, attributed to the emission of the dusty torus. 

For illustration purposes we rebinned each spaxel to 1/3 of its 
original size (0705) to smooth the flux and kinematics maps. We 
applied the same procedure to the continuum image of the standard 
star and noticed that the effect of this procedure in the measured 
FWHM is smaller than 5% and thus the rebinning of the maps does 
not degrade the spatial resolution of the data significantly. 


3 RESULTS 

The left panel of Figure [T] shows an optical image of NGC 2110 
obtained with the Hubble Space Telescope (HST) Wide Field Plan¬ 
etary Camera 2 (WFPC2) through the filter F606W. The FOV of 
the NIFS observations (3" x 3") is indicated by the gray square. 
In the right panel we present a continuum image of the nuclear re¬ 
gion obtained from the NIFS data cube by calculate the average 
of the fluxes between 2.25 to 2.28 pm a region with no emission 
lines. In the bottom panel we present four spectra, extracted within 
a circular aperture with radius of 0725 at the positions marked in 
the K continuum image. The position “+” corresponds to the lo¬ 
cation of the nucleus, defined as the peak of the continuum emis¬ 
sion. The position A [(x,y)=(-072, -072)] was chosen to represent 
a typical extranuclear spectrum, the position B [with coordinates 
(jt,y)=(070, 077) relative to the nucleus] is located 077 north of the 
nucleus and corresponds to a location where a region og enhanced 
H 2 line emission is observed (as seen at the top panel of Fig. 3). Po¬ 
sition C [(jt,y)=(072, 074)] corresponds to a location of enhanced 


Bry emission at 074 northwest of the nucleus. It can be seen that 
the extra-nuclear spectra are very similar (positions A, B and C), 
while the nuclear spectrum (marked with the ’ + ’ sign) is much red¬ 
der, probably due to emission of the dusty torus. At most positions 
the spectra display several emission lines of H 2 , Bry, as well as the 
CO absorption-band heads around 2.3 pm. These lines are identi¬ 
fied in the spectra and have allowed us to obtain the gas and stellar 
kinematics, as discussed below. 


3.1 Stellar kinematics 


We have used the penalised pixel-fitting (pPXF) method of 
ICappellari & Emsellem {2004 ) to derive the stellar kinematics, by 
fitting the 12 COT 229pm, n COA232pm, 13 COT2.34/mi stellar 
absorptions in the K band following the procedure described in 
iRiffel et al. (2008j) . The library of 60 stellar templates of early-type 
stars JWinge et al 1120091) was used to obtain the stellar line-of-sight 
velocity distribution (LOSVD) at each position. 

The top panel of Fig. [2] presents the velocity (V*) map from 
which we subtracted the systemic velocity, obtained by fitting the 
velocity field as discussed in section l4Tl This map shows a rotation 
pattern with blueshifts to the North-West and redshifts to the South- 
East of the nucleus, with a maximum velocity of « 200 km s -1 and 
a major axis along the North-South direction at PA « 170°. The 
bottom panel displays the stellar velocity dispersion (cr*) map with 
values ranging from 50 to 300 km s -1 . White regions in the V* map 
and black regions in the cr* map correspond to positions where the 
signal-to-noise ratio in the CO bands was not high enough to allow 
good fits. At the nucleus of the galaxy, the couldn’t get good fits 
due to the dilution of the CO absorptions by non-stellar continuum 
emission. The dilution may be caused by emission of hot dust, as 
seen in the nuclear spectrum shown in Fig. \l\ The uncertainties 
were derived from the pPXF routine and are smaller than 30 km s -1 
for both the V* and cr* maps. 


3.2 Emission-line flux distributions 

In order to construct maps for the gas flux distributions, centroid 
velocities and velocity dispersions, we fitted the emission lines of 
H 2 and Bry with Gaussian curves using a modifie d vers ion of the 
emission line profile fitting routine PROFIT felffefeOld) . This an 
automated routine written in IDL language and can be used to fit 
emission lines with Gaussian curves or Gauss-Hermite Series using 
the MPFITFUN routine dMarkwardtl2009h to perform the nonlinear 
least-s quares fit. More details of the routine can be found in Riffel 

( 2010f) . 

The flux distributions for the emission lines H 2 A 2.121 8pm 
and Bry T 2.1661 pm (in units of 10 -17 erg s _1 cm -2 ) are shown, 
respectively, in the top and bottom panels of Fig. [3] The 
H 2 T2.1218//m flux map shows emission over the whole FOV, 
with two regions of higher intensity relative to the surround¬ 
ing regions, one at the nucleus and other at « 077 to the 
North of it. The flux distributions for other 7 emission lines 
of H 2 (H 2 T 2.22344 pm, H 2 T 2.24776 pm, H 2 T 2.40847 pm, 
H 2 T 2.41367 pm, H 2 T 2.42180 pm, H 2 T 2.43697 pm and 
H 2 T 2.45485 pm) are similar to that of H 2 A 2.1218 pm, pre¬ 
senting however a lower signal to noise ratio and for this reason 
they are not shown. 

The Bry flux map shows emission only covering part of the 
FOV being extended along the South-East-North-West direction, 
with the peak flux at « 074 North-West of the nucleus. The un¬ 
certainties in flux are smaller than 20% at most locations for both 
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Figure 1 . The top left panel shows an HST image of NGC 2110 with the filter F606W, with the FOV of the NIFS observations indicated by the central square. 
The top right panel shows an average continuum image (obtained in the spectral range between 2.25 - 2.28/zm) in flux units of 10“ 17 erg s -1 cm -2 A -1 ; the 
central cross indicates the position of the nucleus and the other points indicate the positions from which the spectra were extracted. The bottom panel displays 
the spetcra in arbitrary flux units per wavelength, which were extracted within a circular aperture with radius of 0725 and corrected for Doppler shift. 


emission lines and black regions in Figure [3] correspond to loca¬ 
tions where the uncertainties are larger than 50% and no reliable 
fits were obtained. 

Tab. [T] displays the fluxes for the detected emission lines at 
positions A (next to the nucleus), B (the northern blob of enhanced 
emission in the H 2 flux map of Fig. [3]) and C (corresponding to the 
peak of the Bry emission at 074 North-West of the nucleus in the 
flux map of Fig. [3}, integrated within a circular aperture of radius 
0725. 


3.3 Centroid velocity and velocity dispersion maps for the 
emitting gas 

The top panel of Fig. 0| shows the centroid velocity map for the 
molecular gas in units of kms -1 . The systemic velocity, as de¬ 
termined from the fit of the stellar velocity field, was subtracted 
and the cross indicates the position of the nucleus. The H 2 veloc¬ 
ity field shows a rotation pattern with blueshifts to the North (thus 
approaching) and redshifts to the South (thus receding), with a ve¬ 
locity amplitude of « 200kms -1 . But, although presenting similar 
rotation to that observed in the stellar velocity field (Fig. 0, the 
H 2 velocity field presents additional components. Two structures 


Table 1 . Emission-line fluxes (units of 10 -16 erg s -1 cm -2 ) for the A, B and 
C regions identified in Fig. [3] integrated within a circular aperture with ra¬ 
dius of 0725. 


Aq (//m) 

ID 

Nuc 

A (0.2" SE) 

B (0.7" N) 

C (0.45" NW) 

2.12183 

H 2 1-0S(1) 

17.3+4.2 

5.55+0.86 

5.95+0.93 

4.14+0.19 

2.16612 

HIBry 


1.12+0.10 

1.33+0.19 

2.46+0.59 

2.22344 

H 2 1-0 S(0) 

- 

1.14+0.11 

2.14+0.41 

1.41+0.18 

2.24776 

H 2 2-1 S(l) 

- 

0.83+0.08 

0.68+0.05 

0.48+0.05 

2.40847 

h 2 1-0Q(1) 

- 

3.23+0.48 

4.52+0.57 

3.34+0.26 

2.41367 

H 2 1-0 Q(2) 

- 

0.84+0.10 

1.54+0.57 

0.90+0.33 

2.42180 

H 2 1-0 Q(3) 

- 

4.24+0.36 

4.45+0.48 

3.38+0.20 

2.43697 

H 2 1-0 Q(4) 

- 

0.79+0.07 

1.34+0.18 

0.67+0.18 

2.45485 

H 2 1-0 Q(5) 

- 

2.86+0.26 

5.04+0.50 

1.93+0.77 


resembling spiral arms are observed: one to the North in blueshift 
and the other to the South in redshift. Closer in, excess blueshifts 
are observed to the West and excess redshifts to the East, at dis¬ 
tances from the nucleus smaller than 075 (73 pc). The dark gray 
curves with an arrow represent the gas flowing to the central region 
along nuclear spirals. 

Figure [5] shows the velocity dispersion (cr) maps for the H 2 
and Bry emitting gas. Both maps have been corrected for the in- 
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Figure 2. Stellar kinematic maps obtained from the pPXF fit. Top: centroid 
velocity map. Bottom: velocity dispersion map. The color bar shows the 
velocity values in units of km s -1 and the central cross marks the position 
of the nucleus. 


Figure 3. Top: H 2 flux distribution. Bottom: Bry flux distribution. The cen¬ 
tral cross marks the position of the nucleus, the color bars show the range of 
flux values for each emission line in units of 10“ 17 erg s -1 cm -2 and white 
contours correspond to the radio emission. 


strumental broadening. The top panel shows the cr map for the 
H 2 T2.1218jum emission line, with uncertainties of « 10km s -1 . 
High velocity dispersions of « 150 km s -1 are seen in a region ex¬ 
tending along PA^ -45° from North-East to South-West and lower 
values of ~ 60 km s -1 are observed in all other regions. 

The Bry cr map is shown in the bottom panel of Fig. [5] with 
mean uncertainties of 15 kms -1 and presents somewhat higher cr 
values than those observed for the H 2 in the same regions. At most 
locations cr « 100 km s -1 are observed in the regions where the H 2 
cr is ~ 60 km s'. 


3.4 Channel Maps 

Figures [6]and [7]show channel maps along the H 2 A 2.1218yt/m and 
Bry emission line, respectively. In each panel, the flux distributions 
have been integrated within velocity bins of 60 km s -1 , correspond¬ 
ing to two spectral pixels. The flux distributions are shown in loga¬ 
rithmic units and the central velocity of each panel (relative to the 


systemic velocity of the galaxy) is shown in the top-left comer. The 
central cross marks the position of the nucleus. 

Fig. [6] shows that the molecular hydrogen presents emission 
from gas with velocities in the range from -300 to +350 kms" 1 . 
The highest blueshifts are seen at 1" North of the nucleus, and 
at a small region surrounding the nucleus, while the highest red- 
shifts are observed predominantly to the South of the nucleus at 
distances 1" from it and between the nucleus and 072 to the East. 
At lower velocities (from « -100 to +100kms -1 ), the H 2 emis¬ 
sion spreads over a larger part of the FOV with the highest inten¬ 
sity structure moving from North-West to South-East as the central 
velocities change from blueshifts to redshifts. 

The channel maps across the Bry emission line profile are 
shown in Fig. [7] The highest blueshifts and redshifts are observed to 
the North-West and South-East of the nucleus, respectively. There 
is a trend similar to that observed for H 2 for velocities between 
« -100 and + 170 km s -1 : as the velocities change from blueshifts 
to redshifts, the highest intensities move from North-West to South- 
East of the nucleus. 
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Figure 6. Velocity channels along the H 2 emission line in « 60 km s -1 bins (larger than velocity resolution of the data FWHM « 44 km s _1 ) centered on the 
velocities indicated in each panel. The central cross marks the position of the nucleus (peak of continuum emission), the intensities are represented in the color 
scale to the right in logarithmic units and the green line corresponds to the line of nodes. The spiral arms are delineated by the blue curves shown in the H 2 
channel maps indicating the location of the inferred inflows discussed in Sec. |4 a71 




Figure 7. Same as Fig.[6]for Br y. 
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Figure 4. Top panel: velocity field of the H 2 emitting gas. Bottom panel: 
velocity field of the Bry emitting gas, restricted around PA^ -30°/150°. 
The color bars are in units of kms -1 and the mean uncertainty in velocity 
is less than 10kms -1 . The spirals delineated by the gray curves shown in 
the H 2 velocity map mark the location of the inferred inflows discussed in 
Sec. l4~4~Tl In both maps, regions with errors higher than 50% in flux were 
masked out, as for the flux map. 

4 DISCUSSION 
4.1 Stellar kinematics 

[Gonzalez Delgado et alJ ( 2002) presented a study of the stellar 
kinematics of the central region of NGC 2110 based on op¬ 
tical long-slit observations along PA=6° of the Ca n triplet 
(228498,8542,8662) absorption lines. They measured a velocity 
amplitude of « 180kms -1 , with the the line of nodes oriented 
along PA= 163° and velocity dispersion of « 260 + 20kms -1 . 
This_cr value agrees, within the uncertainties, with that obtained by 
iNelson & Whittle! (1 1995h (c r = 220 ± 25 km s’ 1 ) using Mg b lines. 

Ferruit et al.l d2004l) obtained the stellar kinematics from Mg 
and Fe absorptions lines in the optical with the integral field spec¬ 
trograph OASIS of the Canada-France-Hawaii Telescope (CFHT) 
in the spectral range 4760 - 5560A over a 4" x 4" FOV. They found 
a velocity amplitude of « 200kms -1 with the orientation of the 


Figure 5. Top: velocity dispersion map for H 2 with contours of the 3.6- 
cm Very Large Array (VLA) radio map superimposed. Bottom: velocity 
dispersion map for Bry. Both maps are in units of kms -1 . 


line of nodes being consistent with the orientation of the photomet¬ 
ric major axis of the galaxy - PA= 163°. The authors also fitted 
the near-IR CO absorption band heads from long slit spectroscopy 
obtained with the NIRSPEC instrument of the Keck-II Telescope. 
They found a good correspondence between the NIRSPEC and OA¬ 
SIS measurements. 

In order to compare our stellar velocity field with those ob¬ 
tained in the previous studies described above, we have fitted a ro¬ 
tating disk model to the stellar velocity field of Fig. [2]), ass uming 
circular orbits in the plane of the galaxy, as in feertola et all 199 ll) : 

17/ D \t/\ _ t 7 , ARcosCV - x ¥ 0 )sin(i)cos p 6 

V I ) - v /A 

R 2 [sin 2 ('P - %) + cos 2 (i)cos 2 ('¥ - T 0 )] + c 2 fos 2 (if 

where R and T* are the coordinates of each pixel in the plane of the 
sky, V s is the systemic velocity, A is the amplitude of the rotation 
curve at the plane of the galaxy, T'o is the position angle of the line 
of nodes, c 0 is a concentration parameter, i is the disk inclination 
relative to the plane of the sky, p is a model fitting parameter. For 
p - 3/2 the system has a finite mass, as for a Plummer potential 
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and for p = 1 the rotation curve is asymptotically 


The top-panel of Fig. [8] shows the rotating disk model that 
provided the best fit of the observed velocities. In the bottom panel 
we show the residual map, obtained from the subtraction of the 
model from the observed velocities. Although at some locations 
the residuals are larger than the calculated uncertainties, the resid¬ 
ual map does not show any systematic structure, and presents resid¬ 
uals larger than 50 km s _1 only for 13% of the spaxels. We thus con¬ 
clude that the model is a good representation of the stellar velocity 
field. The parameters derived from the fit are: systemic velocity cor¬ 
rected to the heliocentric reference frame V s = 2332 + 15 km s~ l , 
'To = 162° ± 0.6°, A = 306 ± 10.5 km s~\ i = 43° ± 2.5°, 
c 0 = 0.9 ± 0.04arcsec. The location of the kinematical center 
was fixed at X 0 = Y 0 = 0 arcsec and p was limited to range be¬ 
tween 1 and 1.5, which is the range of values expected for galaxies 
feertola et al.ll 199 if) . The best fit was obtained with p - 1.5. The 
systemic velocity and the inclination of the disk are in good agree¬ 
ment with those listed in the NED database. The orientation of the 
line of nodes agrees with the one quoted in Gonzalez Delgado et alii 
,2002). 

We can also compare our results w i th tho se obtained for the 
gas kinematics by ISchnorr-Muller et al] (l2014h using the GMOS 
IFU for the inner l.lxl.6kpc 2 of NGC2110. They identified four 
components in the emitting gas named as: (1) warm gas disk 
(cr = 100 - 220kms _1 ); (2) cold gas disk (<x = 60 - 90km s _1 ); 
(3) nuclear component (cr = 220 - 600 km s -1 ) and (4) north¬ 
ern cloud (cr = 60 — 80kms _1 ). They fitted the velocity field of 
the disk components by a rotation model similar to the one used 
here. For the cold component they found V s = 2309 + lOkms -1 , 
'Fq = 171° + 1° and i = 39° + 1°. For the warm component they 
found V s = 2308 + 11km j - 1 , T'o = 169° + 1°, keeping the in¬ 
clination fixed to the same value of the cold disk. The parameters 
for both disk components are similar to those we obtained for the 
stellar kinematics, as well as the amplitude of the velocity field in 
the inner 175 obtained for the optical emitting gas. 

The average stellar velocity dispersion obtained from the mea¬ 
surements shown in Fig. [2] is cr* = 200 + 20 k ms" 1 , w hich is 
somew hat smaller from those obtain ed bv lGonzalez Delgado et al.1 
d2002h and lNelson & Whittl9 i 19951) using optica l absorption lines 
as mentioned above. However, Rif fel et al.1 J2014h showed that the 
cr* obtained using CO bandheads is systematically smaller than 
the values obtained from the fitting of optical absorption lines 
with a mean fractional difference (< (Tf rac = (cr C o ~ cr opt )/cr opt ) of 
14.3%, with some galaxies showing much higher discrepancies (as 
NGC 4569, with a fractional difference of 40%). Integrating the 
K-band spectra over the whole field of view and using the pPXF 
code, as done in Sec. l3.ll we obtain cr* = 230 + 4km s _1 . Using the 
value of Gonzalez Delga doet al.1 (l2002h we find a fractional differ¬ 
ence of -11%, while using the value of iNelson & Whittlj (l995) 
we get a difference of 4%, T hese differences are of the order of the 
discrepancies found in lRiffel et al.1 i2014f) . 

Using this result in the M. - cr* relation iGraham et all 201 ll) 
we obtain a mass for the SMBH of M. = 2.7 | 3 ' 3 x 10 8 M©, 
which agrees with value found by iMoran et~aD ( 2007 ) using a 
value for the velocity dispersion measured in optical observations. 
T he M« we have found in this work is also consistent with that 
of iDurre & Mould 12014 ) obtained using the kinematics of the 
ionized gas (Hei) to derive the mass enclosed in the inner 56 pc 
(« 4 x 10 8 Mo). 
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Figure 8. Rotating disc model for the stellar velocity field (top) and residual 
map between the observed and modeled velocities. The cross marks the 
position of the nucleus and the color bar shows the range of velocities in 
kms -1 . 


4.2 Gas excitation 

The H 2 molecule can be exc ited by processes such as heat- 
ing of the gas by shock s_ (Hollcnbach & McKee 1 19891) or 
by X-rays ( Maloney. Hollenbach & Tielensl Il996h . known 
as thermal processes, and by non -jher mal processe s such 
as absorption of soft-UV iBlack & van DishoeclJ Il987l) 

and far-UV at clou ds with density larger than 10 4 cm -3 

(ISternberg & Dalgarnol ll989h . The origin of the H 2 emission 
in AGN has been discussed in several recent studies, based 


on long-slit (Store 

li-Bergmann et al. 19991: Reunanen et all 

120021; iDavies et al. 

200j) and inte 

gral field spectroscopy 

(Riffel etaD 2006; Sanchez et al. 2009 

; Storchi-Bergmann et al. 

20091: Mazzalav. Rodriguez-Ardila & Komossal 20lOh. In the K 


band there are several emission lines from H 2 , that can be used 
to study their origin. We used the ratio H 2 T2.2477//12.1218 to 
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assess if the H 2 emission is due to a thermal or non-thermal 
process at A and B locations in Fig. [3] For thermal pro¬ 
cesses, typical values are « 0.1 - 0.2, while for non-thermal 
processes, typi c al val ues are » 0.55 (e.g. M ouri H. (Hi); 

iReunanen et all d2002h : IStorchi-Bergmann et al.l (120091) ). For 
NGC 2110 we find H 2 42.2477/42.1218 « 0.15 for the position A 
and H 2 42.2477/42.1218 ^0.11 for positions B and C, indicating 
that the H 2 is thermalized by shocks and/or by X-ray emission. 
However the signal-to-noise ratio for the H 2 42.2477 emission line 
was not high enough to further investigate the H 2 42.2477/42.1218 
ratio at other locations. 

The lower vibrational energy states of the H 2 molecule can 
also be th ermalized in high gas densiti es even if excited by far-UV 
radiation iSternberg & Dalgamoll 19891) . and long-slit spectroscopy 
of the circumnuclear region of a few active galaxies suggests that 
the this might be the case of some objects dDavies et al .11 20051) . 

Assuming that the H 2 is thermalized, we can find the exci¬ 
tation te mperature by fi tting the observed fluxes by the following 
equation ( IScoville et all 19821) : 


log 



: constant - 


T x 


(2) 


where F x is the flux of the i th H 2 line, A t is its wavelength, A t is the 
spontaneous emission coefficient, g t is the statistical weight of the 
upper level of the transition, T { is the energy of the level expressed 
as a temperature and T exc is the excitation temperature. This equa¬ 
tion assumes a ratio of transitions ortho:para of 3:1 and is valid 
only for thermal equilibrium, what seems to be the case for most 
AGN, as found in the studies discussed above: non-thermal pro¬ 
cesses contribute very little to the H 2 emission, at least in the cases 
that do not show intense star formation. In the case of NGC 2110, 
the H 2 42.2477/42.1218 ratio discussed above supports that the H 2 
is thermalized. 

In Fig. [9] we plot the A upp = F x Ai/A x g x (plus an arbitrary con¬ 
stant) versus E upp = T t for the positions A, B and C shown in Fig.0 
in which ortho transitions are represented with open symbols and 
the best linear fit of equation [2] is shown as a continuous line. This 
result confirms the thermal excitation. Although the uncertainty in 
the fit might be high, since there is only one point at the high en¬ 
ergy region, it suggests that the vibrational and rotational states of 
the H 2 molecule are in thermal equilibrium and this favors thermal 
processes as the main excitation mechanism of the H 2 emission. 
The derived excitation temperatures for positions A, B and C are 
T a = 2644 + 249 K, T B = 2133 + 118 K and T c = 2143 + 140 K, 
respectively. 

In order to further investigate the origin of the 
H 2 emission, we estimate the vibrational tempera¬ 
ture (T vib ) using the H 2 42.1218/2.2477 line ratio in 
7Vj£=5600/ln( 1.355 xFH 2 /12. 1218/F H 2 .12.2477) ( Reunanen et a l 
2002 ) using the fluxes from Tabled We find that T vib is in the range 
2260-2540 K, which are very similar to the values found above for 
the excitation temperature. The fact that the excitation temperature, 
obtained using both rotational and vibrational transitions, and the 
vibrational temperature are similar, supports the assumption above 
that both rotational and vibrational states are thermalized. 

We also investigated the excitation mechanism of the 
H 2 using the H 2 42.1218 /Bry line ratio. In starbursts, 
where the main excitation mechanism is UV radiation, 
H 2 42.1218 /Bry < 0.6, while for AGN, where the H 2 is ex¬ 
cited by shocks, by X-rays from the nucleus or by UV radiation, 
0.6 < H 2 42.1218 /Bry < 2, as found from lon g-sli t spec - 
troscopy for the nuclear aperture feodriguez-Ardila et alJ 120041 : 



Figure 9. Relation between V upp = F x A\ /A\ g\ and F upp = T t for the H 2 
emission lines for thermal excitation at locations A, B and C (from top to 
bottom) in Fig. [3] Ortho transitions are shown as filled circles and para 
transitions as open circles. 


Rodriguez-Ardila, Rilfel & Pastorizal l2005l) . Recent studies using 
long-slit spectroscopy and IFU spectroscopy indicate a broader 
range of H 2 /Bry ratio v alues for AGN, with values i n the range 
0.6-6 dRiffel et al .1120131) and 0.5-8 JColina et al .1120141) . Figure ITOl 
shows the H 2 /Bry line ratio map for NGC 2110, where the values 
are between 0.5 and 8 at most locations indicating excitation by 
the AGN. Approximately in 80% of the points the values are larger 
than 2, but 91% of the points show values between 0.5 and 8, thus 
still typical of Seyfert galaxies. 

Our Bry flux map (Fig. 0 pres ents a fairly similar st r ucture 
to that seen in the [Nn] flux map from lSchnorr -Muller et al. (2014) 
and to that found for the [Oiii] emission by Ferruit et HI J2004), 
where the ionized gas emission is stronger to the North of the 
nucleus than to the South, although the Bry emission is less ex¬ 
tended than both [Nn] and [Oiii]. The Bry emission is collimated 
along the direction South-East - North-West (PA « -30°/150°) 
that is the same o rientation of the strongest ga s emission in the op¬ 
tical as shown in ISchnorr-Mtiller et al.l (120141) . The Bry emission 
can thus be attributed to ionization by the collimated AGN radi¬ 
ation, that illuminates the g alaxy disk in that direction. While in 
ISchnorr-Mtiller et all d2014r) they have detected (much fainter) gas 
emission beyond this region of strongest emission, we are not able 
to detect this emission in Bry, as this line is much less intense than 
the optical lines. 


4.3 Mass of the molecular and ionized gas 

The mass of hot H 2 in the inner 3" x 3" can be estimated as follows 
dScoville et al.ll 19821) : 

2 trip ^h 2 , 12.1218 (F 
™h 2 - — 7 -:—;- 

/v=U=3 As(l)«V 


* 5.0776 x 10 13 


Fh 2 A2.121S 
ergs -1 cm -2 



[M 0 ], 


(3) 


where m p is the proton mass, F H2 42.1218 is the flux of 
H 2 4 2.1218//m emission line, d is the galaxy distance, h is the 
Planck’s constant and v is the frequency of the H 2 line. We 
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Figure 10. Line ratio map H 2 42.1218/Bry 


have _a_ssumed a typic al vibrational temperature of T V ib = 2000 K 
IStorchi-Bergmann et alJl2009h , which implies a population frac¬ 
tion / v =i,j=3 = 1. 22 x 10 -2 and a transition probabilit y A$(i) = 
3.47 x 10 -7 s -1 (iTurner et al.1 1 19771 : iRiffel et al.l l2008h . Integrat¬ 
ing over the whole field of NIFS, the flux is Fu 2 A 2 .ms - 3.1 x 
10“ 14 erg s -1 cm -2 , resulting in a hot molecular hydrogen mass 
of Mh 2 ~ 1.4 X 10 3 M o . This val ue is in good agreement 

with those found for other AGN (Storchi-Bergmann et al. 2009; 
IRiffel. Storchi-Bergmann & Nasaii l20?0al : IRiffel et al. 20131) . The 
amount of cold molecular gas may be much higher than the value 
derived above. iDale et all (2005) found a ratio between the mass 
cold (M co id) and hot molecular in t he rang e 10 -10 7 for a sam¬ 
ple of six nearby galaxies, while Sanchez et al] ( 20091) found a 
ratio 1 x 10 6 for NGC 1068 and iMazzalav et al.l d2013h obtained 
4^coid/4^H 2 ranging from 10 5 to 10 8 . We use her e the conversion 


factor obtained by the most recent study of IMazzalav et al. 1 (1 20131) 
because it includes more galaxies than the previous studies, that can 
be obtained from the expression: 


M co \d 

"mT 


1174 x 


■^H 2 /12.1218 


(4) 


where L H2 A 2.1218 is the luminosity of the H 2 line. We find M co \ d « 
9.92 x 10 8 Mq, that is about 6 orders of magnitude larger than the 
mass of hot molecular gas. 

The mass of ionized gas can be estimated from the measured 
flux of the Bry emission l ine using the following expr ession (de¬ 
rived from the equations in lOsterbrock & Ferland 


Mhh * 3 x 10 17 


Bry 


ergs -1 cm - 


(sk) ! 1M ° l 


(5) 


where F Bry is the integrated flux for the Bry emission line and 
d (30.2 Mpc) is the distance to the galaxy. We assume an elec¬ 
tron temperature T - 10 4 K and electron density N e = 500 cm -3 . 
Integrating over the whole IFU field we obtain F Bry = 6.5 x 


10 1 5 erg cm 2 s \ resulting in M\ 


1.77 x 10 6 M o . 

iDurre & Mould I (820141) . through adaptive optics near-infrared 
integral field spectroscopy with the Keck OSIRIS instrument found 
a total mass of ionized gas of 7.5 x 10 5 M© within a FOV of « 
2.3 x 1.8 arcsec. This value is about 40% of the value we have 


obtained, but this is consistent with the fact that the area covered 
by our FOV is 2.2 times larger. 

The mass of ionized gas derived for NGC 2110 is thus 3 or¬ 
ders of magnitude smaller than the estimated total mass of molec¬ 
ular gas , a result that is i n agreement with those of p r evious 
studies (jFerruit et aL 20041 : IStorchi-Bergmann et 3B 2010 : 
Riffel et al. 2008, 201J - 


4.4 Kinematics of the emitting gas 

The velocity fields of the H 2 and Bry emitting gases (Fig. [4j are 
similar to the stellar one (Fig. O, with blueshifts to the North-West 
and redshifts to the South-East of the nucleus. However, the H 2 
kinematics shows deviations from rotation in regions next to the 
nucleus as well as in regions « 170 to the North and South, that 
show a spiral pattern. 

Although the inclination of the disk and orientation of the line 
of nodes obtained from the fit of the stellar velocity field are no t 
that different from those obtained by ISchnorr-Miiller et all d2014l) 
for the gas in the optical, there are significant differences between 
the H 2 velocity field and the stellar velocity field. Several works in¬ 
deed show that the gas and stellar kinematics are not al ways con sis¬ 
tent with each other in the central region of galaxies (iP avies et all 
2014; Westobv et ah 20121 : Rifle IStorchi-Bergmann, Rifle] 120151 ; 


Stoklasova et al.l 2009th 


In order to isolate the non-circular components of the molecu¬ 
lar gas velocity field, we fitted the H 2 velocity field using the Eq.Q] 
but keeping fixed the geometric parameters of the disk and the sys¬ 
temic velocity derived from the fitting of the stellar velocity field. 
This procedure allows us to derive the maximal rotation, that is ex¬ 
pected to be larger for the gas than for the stars, since the gas may is 
usually located in a thin disk (while the stars have a higher velocity 
dispersion). The resulting best model is very similar to that for the 
stars (Fig. I u , but with a velocity amplitude of A = 473 + 20 km s 1 . 
We have built a residual map between the H 2 velocity field (Fig. 0 
and the fitted best disk model. This H 2 residual map is shown in 
Fig. QT] Assuming that the spiral arms are trailing, we conclude 
that the far side of the galaxy is the North-East and the near side is 
the South-West. 

Beyond the inner 075 (73 pc) of the H 2 residual map we find 
what seems to be a spiral arm in blueshift to the North and North- 
East of the nucleus in the far side of the galaxy, that could be due to 
gas flowing in towards the nucleus. In the near side of the galaxy we 
find mostly redshifts, that can also be interpreted as due to inflows. 

We thus interpret most of the residuals observed in Fig. \TT\ 
as being originated from gas in the plane of the galaxy flowing 
towards the nucleus. This interpretation is supported by the channel 
maps in H 2 shown in Fig. [6] in the channel maps from -245 km s -1 
down to -65 km s -1 , the blueshifted arm can be observed in the far 
side of the galaxy plane, while in the channel maps from 116 km s -1 
to 356 km s -1 , the redshifted arm can be observed in the near side 
of the galaxy plane, supporting inflows towards the line of nodes 
where a structure resembling a nuclear bar is observed in the lower 
velocity channels (blueshifts and redshifts). This bar runs from the 
North-West to the South-East, and is actually a bit tilted relative to 
the line of nodes. 

Within the inner 075 (73 pc), the residual map of Fig. fill 
shows blueshifts in the near side of the galaxy plane and red¬ 
shifts in the far side along PA ^ 90° that we attribute to an 
H 2 nuclear outflow. A compact outflow, but for the ionized gas, 
has been r eported in pr evious optic al long-slit observations by 
iGonzalez Delgado et al.l d2002h . iRosario et al.l (l2010l) . as well as 
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Figure 11. Velocity residual map for H2/12.1218// m obtained from the dif¬ 
ference between the H 2 velocity field (Fig. |4j and the gas velocity model. 
The map is in units of km s _1 . 



by our group in lSchnorr-Muller et al.l d2014h using optical integral 
field spectroscopy. 

Schnorr-Miille r et al.l d2014h has found two kinematic com¬ 
ponents in disc rotation in NGC 2110, that they have called the 
“warm” (higher velocity dispersion) and “cold” (lower velocity 
dispersion) disks (their Figs. 11 and 12). Their residual velocity 
maps for the warm component shows positive values to the North- 
East of the nucleus and negative values to the South-West of it, 
in a very compact structure (within the inner arcsec). They inter¬ 
preted this component as a nuclear outflow and we conclude it 
could be identified with the compact outflow we see in H 2 (Fig.fTTT). 
Similarly, the veloci t y resid ual map for the cold component of 
ISchnorr-Muller et all d2014l) shows some excess redshifts to the 
South-West and blueshifts to the North-East that are similar to the 
residuals found in the H2 kinematics, that we thus identify with the 
inflows we have found in H 2 . 


4.4.1 Feeding of the AGN 


Nu 2 


Mh 2 

2 m v n r 2 d h 


(7) 


We find A^h 2 ~ 6.22 x 10 -3 cm -3 , which results in an inflow rate of 
« 4.55 x 10 -4 M© yr -1 . The inflow rate deriver above is small, and 
might be only a fraction of the total molecular gas flowing towards 
the center of the galaxy, as the warm molecular gas we observe in 
the near-IR may only be the heated skin of a much larger reser¬ 
voir of cold molecular gas, as we have discussed in Sec. 14.31 Using 
the conversion factor for the cold/warm molecular gas obtained in 
Sec. 14.31 we find that the molecular gas mass inflow rate could be 
higher than 10 M© yr -1 , 


4.4 .2 Feedback of the AGN 

As we have also found outflows within the inner 075 (73 pc), we 
can also calculate their power. These outflows have been identified 
as the blueshifts and redshifts in the residual velocity map of the 
molecular hydrogen (Fig. HD along PA « 120° within the inner 
075. 

The outflow observed in H 2 has a similar orientation to that of 
the co mpact outflow seen in ionized gas by ISchnorr-Muller et aD 
(120141) . who also concluded that the outflow is not oriented along 
the ionization cone. Outlfows not co-spatial with the ionization 
cone have been observed for other active galaxies. 

iRiffel et akl d2014bh report an equatorial outflow detected us¬ 
ing NIFS observations for the Seyfe rt 2 NGC 5929 (see also 
IRiffel. Storchi-Beremann. Ritfc] i2015l) . For NGC 2110 we pro¬ 
posed that the outflow is a result from gas in the plane of the galaxy 
being pushed by clouds o utflowing from the nuc l eus. T he same hy¬ 
pothesis was put forth bv ISchnorr-Miiller et al.l d2014h . where the 
authors observed a compact outflow from the nucleus of NGC2110, 
suggesting a spherical geometry, using optical integral field spec¬ 
troscopy. In the case of H 2 , most of the gas is concentrated at the 
disk of the galaxy, and thus the interaction of the expanding spheri¬ 
cal bubble is seen only at locations where the H 2 is present, forming 
the bipolar structure seen in the residual map. 

On the basis of our and other’s previous studies of the mor¬ 
phology of the outflow (as its extent is of the order of the PSF 
FWHM), we assume a bi-conical structure with an opening angle of 
60° with a radius r « 31 pc (consistent with the constraints posed 
by the flux distribution of the outflowing gas) and use the following 
expression: 


Once we have concluded that the molecular gas kinematics suggest 
the presence of inflows towards the center, we can use the velocity 
field maps to estimate the mass inflow rate as follows: 

M H2 =2m f ,N H2 v7rr 2 , (6) 

where m p is the proton mass, v = v 0 b S /seni is the velocity of 
the gas flowing towards the center in the plane of the galaxy, 
v 0 bs = 80 km s -1 is the observed velocity in H 2 the residual map 
and i =_42Ms_the_disk inclination relative to the plane of the sky 
dGonzalez Delgado et al.lll2002 ). This value of i is very similar to 
that obtained by fitting the stellar kinematics (i « 43°). We adopt 
r - 0.3"(~ 44 pc) as the radius of a circular cross section in the 
nuclear arm to the North/North-East of the nucleus at 1" from the 
nucleus. Nu 2 is the molecular gas density and ms = 2 is the num¬ 
ber of spiral arms. We have assumed that the molecular gas is lo¬ 
cated in a disk with radius r d = 1.5" (220 pc) corresponding to the 
NIFS FOV and a typical thickness of h- 30 pc dFticks et al120091) . 

The value of the gas density was obtained from: 


2 m D A H o v n r 2 n 

M out = P H * -. (8) 

sen# 

where v is the outflowing gas velocity, n - 2 accounts for the 
two sides of the bi-cone and 6 is the angle that the central axis 
of the cone makes with the plane of the sky. The velocity can 
be obtained directly from Figure fill and is v = 70 kms -1 . Using 
the estimated 0 = 18° bv [Rosario et al.l (2010), we obtain M out = 
4.3 x 10 -4 M© yr -1 . But again, as for the mass inflow rate, the total 
mass outflow rate, that is probably dominated by the cold molecu¬ 
lar gas, can be 5-7 orders of magni tudes higher. The mass outfl ow 
rate obtained for the ionized gas bv ISchnorr-Muller et al.1 d2014l) . is 
0.9 M© yr -1 , thus much higher than that in the warm molecular gas, 
although it may be comparable to that in the cold molecular gas 
in NGC 2110, as discussed above. This result for NGC 2110 differs 
from those we have obtained for most galaxies we have studied 
so far, where the H 2 gas is usually rotating in the galaxy plane or 
in inflow. It is seldom found in outfl ow. Outflows are more fre- 
quently found in the ionized gas (e.g. ISchnorr-Muller et akll2014l : 
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Riff el & Storchi-Bergmannl 1201 lal : IStorchi-Bergmann et al.1 l20ld : 
Riff el et al .1120131) . 

iDurre & Mould ( 2014 ) report the presence of a nuclear bar 
seen in the [Fe n]flux map, but the velocity field observed in [Fe n] 
shows only blueshifts, although this depends on the adopted sys¬ 
temic velocity. The position angle of this bar, as well as its extent 
(~90 pc) is consistent with what we have interpreted as an outflow 
seen in our H 2 velocity map, that shows blueshifts to one side of the 
nucleus but redshifts to the other side. If this structure is a bar, one 
possibility is that the inflows are feeding the bar, whose orientation 
(P.A.~ 125°) is close to the galaxy major axis (P.A.~ 162°). Unfor¬ 
tunately we do not have observations in the J or H band to compare 
the [Fen]emitting gas kinematics with that of the H 2 emitting gas. 
Further observations and modeling would be necessary in order to 
settle on the nature of this feature. 


4.4.3 AGN accretion rate 

The inflow and outflow rates estimated above can be compared to 
the accretion rate onto the SMBH that is needed to power the AGN. 
The mass accretion rate to the AGN can be estimated via: 


C 2 T] ’ 


(9) 


where c is the speed of light, rj is the efficiency of conversion of 
matter into energy (for Seyfert galaxies tj « 0.1), and L bo i is the 
AGN bolometric luminosity. The X-ray luminosit y of the AGN in 
NGC 2110 is estimated as L x = 2.9 x 10 42 erg s -1 ([Pellegrini et al.l 
20101) , and in order to obtain the bolometric luminosity we adopt 
the usual relation L bo i = 10L*. The resulting accretion rate is thus: 
m- 1.6 x 10 _3 M©yr _1 . 

_If we use instead the bolometric luminosity from lRosario et al.l 

(hold ), of L bol = 5 x 10 43 erg s’ 1 , we obtain a somewhat larger value 
of m - 2.8 x 10 -3 M 0 yr -1 than that obtained above. 

The accretion rate is approximately one order of magnitude 
higher than the inflow and outflow rates in warm molecular gas, 
indicating that it is necessary more mass flowing towards the center 
to feed the SMBH, but, as pointed above, the inflowing mass is 
probably dominated by cold molecular gas. 

The mass inflow rate in warm H 2 obtained for 
NGC 2110 is_ similar to those obtained for other Seyfert 
galaxies Jstorchi-Bergmann jt aD l20ld: iRiffel et al.l 12004 
IRiffel. Storchi-Bergmann & Naaail 2010ah with similar luminosi- 
ties to that of NGC 2110. The outflow rate in warm molecular gas 
is about three orders of magnitude smaller than mass outflow rates 
we have observed in ionized gas for other Seyfert galaxies (e.g. 
IRiffel et al1l2008L 120131 : ISanchez et al .1120091) and up to 6 orders of 
magnitude smaller than those observed for high luminosity Seyfert 
2 galaxies (e.g. iMcElrov et al.l 120141) . But the compact outflow 
we have observed in warm H 2 for NGC 2110 may represent only 
a small fraction of the total outflow from the nucleus of this 
galaxy, that is probably also dominated by cold molecular gas. 
A possible scenario for the inner region of NGC 2110 is that the 
inflows observed in H 2 feed the central region of the galaxy and 
the accreted gas then triggers both star formation and the nuclear 
activity. Indeed, young stellar cluster s have been observed around 
the nucleus of NGC2110 bv lDurre & Mould I (120141) . 


5 CONCLUSIONS 

Using integral field spectroscopy in the near-infrared obtained 
with the Gemini instrument NIFS, we have mapped the stellar 
and gaseous kinematics within the inner ^250 pc of the Seyfert 
2 galaxy NGC 2110, as well as the ionized and molecular gas flux 
distributions at an angular resolution 0.15". Our main conclusions 
are: 

• The stellar kinematics shows velocity dispersion reaching 
« 250 km s~ x around the nucleus as well as a rotation pattern with 
a similar velocity amplitude; 

• From the relation M. — cr*, we estimate the mass of the SMBH 
as M. = 2.7 ! 3 ; 3 x 10 8 M o . 

• The flux distributions of the H 2 and Bry emission lines are 
different: while the Bry flux seems to be collimated along PA « 
-30° (from South-East to North-West), the H 2 emission is observed 
over the whole FOV: 

• From the emission-line ratios, we conclude that the H 2 emis¬ 
sion is due to thermal processes (heating X-rays from the AGN 
and/or shocks), with an excitation temperature in the range « 2100- 
2700 K. 

• The warm molecular gas and ionized gas masses within the 
observed FOV were estimated as M H2 ~ 1.4 x 10 3 M© and M H n » 
1.7 x 10 6 M 0 , respectively; 

• The velocity field of the gas shows a rotation pattern similar 
to that observed for the stars. However, deviations from this pattern 
are observed, due to both inflows and outflows, more clearly ob¬ 
served in the molecular gas emission, as the ionized gas emission 
covers just a small part of the FOV; 

• Inflows in warm H 2 along spiral arms are observed beyond the 
inner « 70 pc, at a mass inflow rate of 4.6 x 10 -4 M© yr -1 ; 

• Outflows in warm H 2 are observed within the inner « 70 pc, at 
a mass outflow rate of 4.3 x 10 -4 M© yr -1 . 

The uncertainties in the above mass flow rates might be high 
due to assumptions about the gas density and geometry. They are 
one order of magnitude lower than the accretion rate to the AGN. 
But we point out that these rates are in warm molecular gas, that 
may be only the heated skin of a much larger cold molecular gas 
reservoir. If we consider that previous observations of both warm 
and cold molecular gas in nearby AGNs reveal that typical ratios 
between the cold and warm molecular gas masses are in the range 
10 5 - 10 7 , the above flow rates may be as high as a few to tens of 
solar masses per year. Similarly, the total mass of molecular gas 
may reach 10 9 M©. 
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